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Abstract 



^ ■ The paper deals with a theoretical study of the transport of a fluid in a channel, which 

takes place by the phenomenon of peristalsis. A mathematical analysis of the said prob- 



lem has been presented. The analysis involves the application of a suitable perturbation 



l> 
in 

j->^ . technique. The velocity profile and the critical pressure for the occurrence of reflux are 

investigated with particular emphasis by using appropriate numerical methods. The effects 
of various parameters, such as Reynolds number, pressure gradient, porosity parameter, 
Darcy number, slip parameter, amplitude ratio and wave number on velocity and critical 
pressure for reflux are investigated in detail. The computed results are compared with 
a previous analytical work and an experimental investigation reported earlier in existing 
scientiflc literatures. The results of the present study are in conformity to both of them. 
The study has got some relevance to the physiological flow of bile in the common bile duct 
in a pathological state. It reveals that in the presence of gallstones, bile velocity increases 
as the value of the porosity parameter increases, while the critical pressure for reflux de- 
creases as porosity increases. 
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1 Introduction 

Peristaltic pumping [1, 2] of physiological fluids induced by a progressive wave of area contraction 
or expansion along the length of a distensible tube has drawn serious attention of researchers 
working in the area of physiological fluid dynamics. Movement of various physiological fluids, 
such as the transport of urine from the kidney to the bladder through ureters, movement of food 
material through the digestive tract, transport 
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a 


Traveling wave amplitude 


C20,C21,C22 


Integration constants 


d 


Half width of mean spacing between solid boundaries 


D 


Constant defined in equation (34) 


m 


Function defined in equation (35) 


k 


Permeability parameter/Darcy number 


p 


Fluid pressure 


dp 
dx 


Pressure gradient averaged over a period of time 


R 


Reynolds number of the fluid 


s 


Slip parameter at the boundary 
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Time 


u,v 


Velocity components in X and Y directions respectively 


x,y 


Rectangular Cartesian co-ordinates 
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Wave number 
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Complex number defined in equation (31) 
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Porosity parameter 
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Wave length of the travelling wave motion of the wall 
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Dynamic viscosity of the fluid 


Hi 


Apparent viscosity of the fluid 





Phase difference 
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Density of the fluid 


01,020,022 


Functions defining ^i,'02 in equation (20)and (21) respectively 
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Stream function 
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Amplitude ratio 


r] 


Vertical displacements of the wall 



of semen in the vas differences, fluids movement of lymphatic fluids in lymphatic vessels, flow 
of bile (cf. Daniel et al. [3]) from gall bladder into the duodenum, spermatozoa in the ductus 
efferentes of the male reproductive tract and cervical canal, movement of ovum in the fallopian 
tube as well as the cilia movement and circulation of blood in small blood vessels, — are all 
performed by the mechanism of peristalsis. The propulsion of some industrial fluids are take 
place by this mechanism. Peristaltic transport has also found various applications in roller and 
flnger pumps, heart-lung machines, blood pump machines, dialysis machines and also transport 
of noxious fluid in nuclear industries. 

A few theoretical studies on peristaltic transport of different physiological fluids were carried 
out by Usha and Rao [4], Mishra and Rao [5, 6], Misra et al [7, 8, 9, 10, 11, 12], as well as Eytan et 
al. [13]. The analyses were mostly restricted to consideration of small peristaltic wave amplitude 
and the assumption that the fluid inertia is negligible. 

Taylor [14] conducted an investigation on asymmetric wave propagation in wavy sheets 
with the main objective of deriving some information regarding the mechanical interaction be- 
tween spermatozoa. In order to study some fluid dynamical aspect of the problem of peristaltic 
transport in an asymmetric channel under Stokes flow conditions, the boundary integral method 
along with a suitable numerical technique was employed by Pozrikidis [15]. Using the assump- 
tions of thin shell and lubrication theories, Carew and Pedley [16] put forward a mathematical 
analysis for the development of peristaltic pumping in the ureter by using similar assumptions. 
By taking into account the wall deformation of the pipe, Antanovskii and Ramkissoon [17] also 
used lubrication theory in order to analyse the peristaltic motion of a compressible viscous fluid 
through a pipe for situations where the pressure drop changes with time. 

Bile (alternatively called as gall) is a greenish yellow secretion that contains various biochem- 
ical substances like bile acids, bile salts, pigments, cholesterol, phospholipids and electrolytic 
chemicals. Bile is produced in the liver. In adult humans, the quantity produced in a day is 
about one litre. After passing through several bile ducts, which penetrate the liver, bile flows 
into to the the common bile duct. It helps accelerate the fat absorption process and there by 
plays an important role in absorbing the vitamins D, E, K and A that are soluble in fat. 

It has been suggested by recent physiological researchers (cf. Vries et al.[18] ) that uter- 
ine peristalsis resulting from myometrial contraction can take place in both symmetric and 
asymmetric directions. Various investigators have carried out different studies pertaining to 
the gastrointestinal tract, intra-pleural membranes, capillary walls, human lung, bile duct, gall 
bladder with stones and small blood vessels, as well as flow in porous tubes and deformable 



porous layers. Keener and Sneyd [19] reported that gastrointestinal tract is surrounded by a 
number of heavily innervated muscle layers which are smooth muscles consisting of many folds. 
There exist pores in the junctions between them, although the junctions are tight. As mentioned 
by Bergel [20], the capillary walls are surrounded by flattened endothelial cell layers, which are 
porous. Li et al. [21] reported that an impulsive magnetic field can be used as a theraptic means 
to treat patients who have stone fragments in their urinary tract. 

Functions of the human biliary system that consists of an organ and a ductal system are to 
create, transport, store as well as release bile into the duodenum to assist digestion of fats. This 
system contains the liver, gallbladder and biliary tract namely cystic, hepatic and common bile 
ducts. Although several analytical and physiological aspects of the human biliary system have 
been analyzed elaborately, we have only insufficient information about the mechanism of bile 
flow in the system. Torosoli and Ramorino [22] reported that pressures in the biliary tree vary 
from 0-14 cm H2O (1 cm H2O = 100 Pa) in the resting gall bladder to approximately 12-20 cm 
H2O in the common bile duct. 

Cholelithiasis is a disease that concerns formation of Gallstones. It has become a major 
health problem worldwide, particularly in adult population. Incidence of the gallstone disease 
indicates considerable geographical and regional variations [23, 24, 25, 26, 27, 28]. 

Lauga and Stone [29] made an experimental attempt to evaluate the effective slip length of 
the resulting flow as a function of the degrees of freedom describing the surface heterogeneities, 
namely the relative width of the no-slip and no-shear stress regions and their distribution along 
the pipe. They gave a possible interpretation of the experimental results which is consistent 
with a large number of distributed slip domains such as nano-size and micron-size nearly flat 
bubbles coating the solid surface. In addition they also suggested that the possibility of a 
shear-dependent effective slip length. 

There are, however, only a few studies on the rheological properties of human bile. Gottschack 
and Lochner [30] examined 33 samples and reported that post-operative T-tube human bile is 
a Maxwell fluid. Coene et al. [31] claimed that 11 of 36 hepatic bile samples displayed non- 
Newtonian behaviour. Lou et al. [32] suggested on the basis of their preliminary measurements 
of fresh human bile collected at 37°C for a healthy person without gallstones that bile is a 
Newtonian fluid. 

From the physiological fluid dynamics point of view, we present here a theoretical study that 
has relevance to the problem of bile transport in the common bile duct in the presence of stones. 
Peristaltic transport of the fluid in a porous channel has been investigated. The symmetric 



sinusoidal peristaltic wave train on the channel walls, impermeable boundary and slip boundary 
condition of Saffman type have been considered in this study. The overall aim of the study has 
been to examine the role of fluid dynamics in the human biliary system, in particular for flow 
in a common bile duct with/without stones. It is expected that the results presented here will 
serve as fairly good theoretical estimates of various prospective fluid mechanical flow governing 
parameters related to the peristaltic transport of bile. 

2 Problem Formulation 





Fig. 1: Geometry of the problem 
In a pathological state, numerous stones are formed in the bile (physiological fluid). In such a 
situation the mixture of the fluid (bile) and the solid particles (stones) forms a dense porous 
mass. We are interested here to study the motion of such a fluid (dense mixture), by considering 
it as a porous mass and treating it as an incompressible viscous fluid. Slip boundary conditions of 
Saffman's type have been considered in the study. We use Cartesian coordinates (x,y), where x is 
measured in the direction of wave propagation and y in the direction normal to the mean position 
of the corresponding common bile duct (cf. Fig.l). Let u, v denote the velocity components in 
the directions of increasing x and y respectively; p, p, v and t stand for the density, pressure, 
kinematic viscosity of the fluid and time respectively. A sinusoidal wave train is propagating 
with a constant speed c along the channel wall, such that r] = acos{^{x — ct))where y = d + rj 
and y = —d — rj represent the upper and lower boundaries of the channel and a is the wave 
amplitude, A the wave length, 2d the width of the channel. 

In order to study the motion of the fluid in the porous structure, we make use of Brinkman's 



equations. These equations and the equation of continuity may be put as 

fdu du du\ dp I d'^u d'^u\ jj, , . 

\dt dx dy J dx ydx"^ dy'^ j k 

I dv dv dv\ dp I d'^v d'^v\ fx , , 

\dt dx dy J dy \dx'^ dy"^ J k 

where /xi = - and e, k are respectively the porosity, permeabihty parameters. 
Introducing the stream function ip, we write 

« = f. - = -!!• 

Ehminating p, one can write the equation that governs the flow of the fluid in terms of ip in the 
following form : 

Ytyy ~r Ytxx ' Yy'ryyx WxWyyy ' YyWxxx WxYxxy \ryyyy ' '^Wxxyy ~r yxxxx) 

~ T^yy ~ T^xx, in which v = - (4) 
k k p 

In the sequel, we shall use the following non-dimensional variables deflned by 

X _ y _ u _ V a _ p - ct _ T] 

X =-7, y = ^, u=-, v = -, e=- P=^, ^ = T7' ^=1' 
a a c c a pc^ a a 

cd - if) - k 2Tid 

i? =— , ^ = — , A:=— , a = — - (5) 

V CCL d A 

If we now drop bars over the symbols, equation (4) may be expressed as 

^yWtyy ~r Ytxx ' YyWyyx WxWyyy ~r WyWxxx WxWxxy ) ~\Wyyyy ' ^Wxxyy ' Wxxxx) 

- -^i'yy--^^xx (6) 

2.1 Boundary Conditions 

Let us now consider the symmetric motion of the flexible walls. The boundary conditions for 
the present problem may be stated as follows : 

(i) non zero velocity slip of Saffman type : u = TsUy or ipy = Tsipyy at y = ±d ± r] 

dri ^TTdc { 27r \ 

[it) symmetric motion of the wall : v = ±t^ or — i/j^ = i — - — sin [^{x ~ ct)j at y = ±d ± rj 

(7) 



3 Solution Procedure 

In order to solve the governing equation along with the boundary conditions, let us express 
the stream function ■?/' as a power series in terms of amplitude ratio e (considered as a small 
quantity), as follows : 

'0 = ^o + e^i + eV2 + --- (8) 

Similarly, we write the pressure gradient -^ as 

ox \OX J ^ \OX J ^ \OX J ^ 

In equation (9) the first term on the right corresponds to an imposed pressure gradient, while 
the higher order terms are conceived of as those arising out of the peristaltic mechanism. Substi- 
tuting (8) into (6), equating the coefficients of like powers of e and neglecting cubes and higher 
power of e, we get 

d 11 

d 11 

(12) 
in which V2 = ^ + ^ anrf ^„ = ^, ^P^y = ^, 2 = 0,1,2. 

Substituting of (8) into (7), we obtain the following boundary conditions by equating the coef- 
ficients of like power of e. 

^o(±l) = and V'oy(±l) = Ts^oyy{±l) (13) 

ipiy{±l) ± i)oyy{±l) COS a{x - t) = Ts{ipiyy{±l) ± ifjQyyy (±1) COS a {x - t)} (14) 

ipix{^^) ± "00x2/(^1) cosa(a; — t) = ^sma{x — t) (15) 

V^2y(±l) ± ipiyy{±l) cos a{x - t) + -cos^a(x - t)'il)oyyy{±l) 

= Ts{ll^2yy{±l) ±l{jlyyy{±l)cOSa{x -t) + -COS^ a{x - t)l/joyyy{±l)} (16) 

7/^2x(±l) ± V^ixj/(±l)cosa(x - t) + - cos^ a{x - t)ilJoa:yy{±l) = (17) 



Considering symmetry and a uniform pressure gradient in the x-direction, the solution of the 
first set of differential equations lead to the following classical Poiseuille flow equation that may 
be written as 

iJo = Ao{-Coy + smh^,} (18) 

yjhere Co = ^ cosh ^ + ^ sinh ^, ^o = —0^. (19) 

y/k \/k k ^Jk Co 

Let us take the solutions of the differential equations (11) and (12), which satisfy the respective 
boundary conditions (14), (15) and (16), (17) given by 

2^1 (x, y, i) = 0i(?/)e^"(^-^*) + 0t(y)e-'"(^-^*) (20) 

2M^, y, t) = <p2o{y) + 022(l/)e^^"(^-^*) + 0;2(?/)e-''"(^-^*\ (21) 

where the asterisk sign denotes the complex conjugate of the corresponding quantity. We now 
substitute (20) and (21) into the differential equations (11) and (12) and the boundary conditions 
(14), (15) and (16), (17). Thus we obtain 

with 0;(±1) + ^ sinh ^ = TS(ly'; - ^^ cosh y^, (23) 

(/)i(±l) = ±1 (24) 

and 0^0-^020 = ^ (0i<^r - '^i0i) (25) 



|^--^-| — (--»(-»^|--f))){|. 



+ ^^ cosh ^022 = -Y- ('/'>i - '/'I'^i ) (26) 



along with the conditions 
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</>',o(±i) ± \ (0;'(±i) + 0r (±1)) + ^ cosh ^ 



= Ts |0'2'o(±l) ± I (0r(±l) + 0r(±l)) + ^ sinh ^1 (27) 

022(±1) ± ^ (01 (±1) + ^r (±1)) + ^ cosh ^ 

= Ts |0'2'2(±1) ± I (0r(±l) + 0r(±l)) + 4^ sinh ^1 (28) 

022(±1) ± ^0'i(±l) = (29) 

In the above- written equations, primes denote derivatives of the function with respect to 'y' ^-^d 
R denotes Reynolds number. 

By solving the set of differential equations presented above along with the respective bound- 
ary conditions, it is possible to obtain the solution of the problem (up to the second order in 
e). In the general case, it is difficult to solve the fourth-order differential equations analytically. 
It is, however, possible to find a closed form solution of the equation for the particular case 
of pumping of an initially stagnant fiuid when there is no imposed pressure gradient, that is 
(af ) = 0- I'^ this case, Aq = and the other coefficients are constants. Such a situation occurs 
in practice, when the fluid is stationary and there is no travelling peristaltic wave. Moreover, 
the maximum pressure gradient that can be created by a small amplitude is of order e^. In 
the pumping range , the zeroth order mean pressure gradient must disappear. Therefore, the 
consideration is not much restrictive. 

Considering Aq = 0, the solution of equation (22), subject to the boundary conditions (23) 
and (24) is found as 

01 (y) = Cii sinh ay + C12 sinh /3y , (30) 



in which 



Cii = 
C12 



/? cosh (5 + s/3^ sinh /3 



a cosh a sinh j3 — /3 sinh a cosh /3 + s(a^ — l3^) sinh a sinh /3 

a cosh a + sa"^ sinh a 
a cosh a sinh P — (3 sinh a cosh /3 -|- s{a^ — (3"^) sinh a sinh /3 

where f3'^ = a"^ + - — 2iaeR, (31) 

k 

when —d — r]<y<d + ri. 

In order to determine the mean flow, we would need only the term 02g in the expansion of 
■?/'2 given by (21). Substituting (30) into (25) and integrating w.r.to 'y' w^ obtain 



1/// _ ^jj _ ^(^J^ 

^20 tV'20 ~ 2 

+2iaeRCi2Cl2 sinh l3y sinh l3*y] + 2C20, (32) 



iaeR — — j Cl^Cu sinh ay sinh jSy + i iaeR + — j C11CI2 sinh ay sinh /3*y 
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where C20 is an arbitrary constant. The equation (32) further leads to 



^20 



^^^^^^ {D-f{l)-sm}cosh^ 2C,ok ^^ 



cosh^ 



y/e I s^e ^ 



v^ 



-«h^ + ^--h^ 



cosh 



in which 

D = (/)'2o(±l)±S02o(±l) 






sinh 






(33) 



1 r 



a\Cii + Cii) sinh a + ^'Cu sinh /? + /3*^Ci2 sinh /3* 
+s {a^iCii + C*i) cosh a + /^^'Cis cosh/3 + /3*3c*2 sinh/3*}] (34) 



tt'^'i f{y) 



iaeR 



iaeR 



^^ n* n \ cosh(a + /3)y cosh(a - /?)?/ 
^11^12 



A; 



(a + /?)^ 



a 



py 



e\ fcosh(a + /?*)?/ cosh(a - /3*)?/' 

■cosh(/3 + /?*)?/ cosh(^-/3*)2/l 



+2mei?Ci2Ci. 
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(35) 



J/3 + /3*)2-f (/3-/3*)'-f, 

One may note that the constant C20 which is proportional to the mean pressure gradient remains 

arbitrary in the solution. In the equations (1-2), if each term is averaged over an interval of 

time i.e the period of oscillation, the solution of the present problem can be derived from (8), 

(18), (20), (21), (30) and (33). The mean pressure gradient may then be expressed as 

^2 ^2 ^2 . ^2 ^ 

(36) 



dp 
dx 



dp' 



dx 



2eR 



ho 



kR 



ho + J*" 



h<^l" 



^i0i) 



+ 0(e^) = — C2o + 0(e^ 
eR 



As a natural consequence of (36), we can write 



dp 
dx 



C20 
'eR 



(37) 



It is now clear that C20 is proportional to the time averaged pressure gradient arising out 
of the peristaltic motion of the fluid. Its value can be determined by employing the conditions 
for a given physiological problem. Noting that i^) is independent of the y-coordinate, the 
expression for the mean velocity in the axial direction can be put in the form 



u 



ho 



m 



[D-fil)-sf'{l))cosh^ 2C2ok 



2 "^^ 2 
Using (36), we can alternatively write 



posh v^ + £v^ sinh ^ 



1 - 



cosh^ 



-^h^ + ^--h^. 



^=2" 



fiy) 



{D -m-sf {l))cosh^ ^^Jdp 



-^h^ + ^--h^ 



-2kR\jf- 

\ dx , 



cosh ^^ 



-^h^ + ^--h^. 



(3^ 



(39) 
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Fig. 2: Variation of D with R, a, Darcy number, porosity and slip parameter 
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Fig.?, FigA 

Fig. 3-4: Mean velocity perturbation function F(y), Fig3: Comparison of our result with Fung 
and Yih [2] for R=l, a =0.4, k=10000, e=0.99, s=0.00001; Fig4: Distribution of F(y) for dif- 
ferent values of R with a =0.25, k=1000, e=0.9, s=0.00001 



For a problem such as the one under our present considerations, it is worthwhile to determine 
the critical reflux condition. This is the condition that determines the location of the point on 
the central line (y=0) at which the mean velocity u{y) vanishes. Using (39), the critical reflux 
condition for the problem is derived in the form 



dp 
dx 



1 



2 critical reflux 



2kR 



/(0)Hh^ + ^--h^ 



D-f{l)-sf'{l) 



-^h^ + ^--h^-l 



(40) 



Use of this critical value enables us to study the mean velocity distribution u{y) at the critical 
reflux condition. 

4 Numerical Results and Discussion 

In the last section, we have presented the analytical expressions of the velocity, time-averaged 
velocity and critical pressure for reflux for the problem under consideration. In this section, we 
intend to investigate the problem and present the computational results for the said quantities 
graphically. Computation has been performed by making use of available experimental and 
theoretical data of bile. The different parameters involved in our forgoing theoretical analysis 
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are Reynolds number R, pressure gradient, porosity parameter e, Darcy number k, slip parame- 
ter, amplitude ratio and wave number. Computation of the time-averaged velocity and critical 
pressure for reflux has been carried out by extensive use of the software Mathematica. 



4.1 Average Velocity in Bile Transport 

Numerical computation based on equation (39) reveals that the averaged axial velocity of bile 
(if treated as a Newtonian fluid) flowing in a porous medium is dominated by the constant 

D , „. JdvS (. ^osh^ 



and the parabolic distribution term 2kR 



^°«h^ + ^^-h^ V^^AV -osh^ + ^sinh^^ 



In addition to these two terms, there is an expression 



V^y 



(/(l) + s/'(l))cosh^ 

that represents the perturbation of the velocity across the channel. Its distribution controls the 
direction of peristaltic mean flow across the cross-section of the bile duct. It is evident from (33) 
and (31) that the constant D that gives the value of 02o ^^ ^^^ boundary depends on Reynolds 
number, porosity, Darcy number, slip parameter and wave number. It arises from the expression 
of the radial component of the first-order axial velocity gradient and is connected with the offset 
boundary condition (27). It may be mentioned here that the no-slip/velocity-slip applies to the 
wavy wall and not to the mean position of the wall. Figs. 2(a-e) give the variation of D with 
the parameters Reynolds number, porosity, Darcy number, wave number and velocity-slip. 

cosh^ 
Furthermore, the parabolic mean velocity distribution term — 2kC2o ' ^ 



-^h^ + ^--h^. 



arises out of the time — averaged second order pressure gradient — — = . 

\dx J ^ eR 

We define F(y)= ~-^^G{y), G(y) being given by (41). Fig. 3 depicts the nature of the function 
F(y). With an aim to validate our results, we have presented in the same figure the results 
reported by Fung and Yih [2] who did not account for the Darcy number, porosity and velocity- 
slip. It is needless to mention that the results presented in this figure on the basis of our study 
correspond to the case when Darcy number, porosity and velocity-slip are neglected (for the 
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Figs. 5-8: Mean velocity perturbation function F(y) (Fig5.) for different values of a with 
R=10, k=1000, e=0.9, s=0.00001; (Fig6.) for different values of k with R=10, a=0.25, e=0.9, 
8=0.00001; (Fig7.) for different values of e with R=10, a=0.25, k=1000, s=0.00001; (Fig8.) for 
different values of s with R=10, a=0.25, e=0.9, k=1000 
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Figs. 9: Time-averaged mean axial velocity profiles for different values of R when a=0.25, e=0.9, 



e=0.1, (|2J = -2.5, k=1000, s=0.0001 
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Figs. 10: Time-averaged mean axial velocity profiles for different values of f #) when R=15, 
a=0.25, e=0.9, e=0.1, k=1000, s=0.0001 

sake of meaningful comparison/validation). It may be noted that the results of both the studies 
have a very close resemblance. In Fig. 4 we have presented the nature of variation of F(y) for 
different values of the Reynolds number for a particular set of values of a, k, e and s. This 
figure shows that F(y) decreases as the Reynolds number R increases. This observation is in 
conformity to that reported in [2]. It may be noted from Figs. 5 and 7 that in the ranges of 
values for a (wave number) and s (slip parameter) considered here, the change in the function 
F(y) is not very much significant. 

Figs. 6 and 8 reveal that with the increase in the Darcy number k and the porosity parameter 
e, the perturbation term F(y) increases. The results correspond to the case of bile containing 
stones. Figs. 9-15 illustrate the individual contribution of Reynolds number, pressure gradient, 
amplitude ratio, Darcy number, the porosity parameter, slip parameter and wave number on the 
time-averaged mean axial velocity profile. It is found that the mean axial velocity increases with 
an increase in Reynolds number (See Figs. 9). It is important to note (cf. Figs. 10) that reflux 
occurs in the central region when pressure gradient attains a certain critical value 0.220966. 
Fig. 11 illustrates that the bile flow is strongly dependent on the amplitude ratio of the wave. 
It reveals that bile velocity is enhanced with the increase in amplitude ratio and that flow 
reversal takes place near the boundary for higher values of the amplitude ratio. These results 
are in conformity (cf. Figs. 9.1-4) to those of an experimental study conducted by Sugita et al. 
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[35]. They investigated pseudolesion of the bile duct caused by flow effect and observed that 
at l.Omm/sec. intervals flow speed of bile ranges between 1.0 and 20 mm/sec. approximately. 
It may be noted from Fig. 12 that although in the central region of the channel, the velocity 
reduces with the decrease in the Darcy number, the velocity profiles are quite different from the 
familiar Poiseuille profile. This shows that the bile velocity decreases as the number of stones 
increases. 

Fig. 13 depicts the variation of the velocity profiles with the porosity parameter e for fixed 
values of Reynolds number, pressure gradient, slip parameter, amplitude ratio and wave number. 
The velocity is found to increase, as the porosity parameter increases in the case when the Darcy 
number exceeds the value 0.05, though its parabolic nature changes for small Darcy number. 
Thus in the situation as the number of stones decreases, the bile velocity gradually increases and 
in the absence of any stone, the bile velocity will be the greatest. The occurrence of a reverse 
trend is observed when the Darcy number equals 0.05. Fig. 14 shows that velocity-slip has a 
strong influence on the axial velocity in the porous case. From Fig. 15, it is observed that wave 
number affects the axial velocity more prominently in the vicinity of the boundary. 
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(c) id) 

Figs. 13: Time averaged mean axial velocity profiles for different values of e when R=15, e=0.1, 
a=0.25, (If) = -2.5, s=0.0001 (a) k=l; (b) k=0.5; (c) k=0.1; (d) k=0.05 




(a) (b) 

Figs.: 14: Time averaged mean axial velocity profiles for different values of s when R=15, e=0.1, 



0.25, (|£) = -2.5, e=0.7 (a) k=l; (b) k=0.05 



dxJ2 




Figs. 15: Time averaged mean axial velocity profiles for different values of a when R=15, 



e=0.1, (g) = -2.5, e=0.7, k=100, s=.001 
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4.2 Critical Pressure for Reflux 

It is known that bacteria and some other materials sometimes move from the bladder to the 
kidney or from one kidney to the other in the direction opposite to the direction of urine flow. 
This phenomenon is referred to as 'ureteral reflux' by physiologists. Severity of diseases such 
as tuberculosis, interstitial cystitis, duct stone etc. is often enhanced due to this reflux. It 
may also happen that due to reflux, from the common bile duct, bile flows into the gallbladder 
and is stored there. Results for the critical pressure for reflux R(-ff-) computed 

\ox / 2 critical reflux 

from equation (40) are displayed in Figs. 16-20. Results presented in Figs. 16 match with 
those of Fung and Yih [2]. This figure illustrates that the said critical reflux decreases as the 
Reynolds number increases. It is observed from Figs. 17-19 that the critical pressure reduces 
with an increase in wave number, whereas it first increases with the increase in Darcy number 
and subsequently maintains nearly a constant value. It is also worthwhile to observe that the 
porosity parameter significantly affects the magnitude of the critical pressure. From Fig. 18, 
it can be conjectured that when the number of stones in the bile is very high, reflux occurs 
when the critical pressure is relatively small. The effect of slip parameter on the critical reflux 
is revealed in Figs. 20. 
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Figs. 19: Effect of e on (^) R=20, a=0.2, k=10, s=0.001 

\ox J 2critical reflux 

5 Summary and Conclusion 

The peristaltic transport of a fluid has been investigated here. The flow is considered to take 
place in a porous channel. The study pertains more particularly to a situation where the 
Reynolds number is low and the curvature of the channel is quite small. The speciflc problem 
has got relevance to the physiological problem of flow of bile in the common bile duct (cf. 
[3, 25, 26, 27, 35, 36, 37, 38, 39]). As mentioned in Section 1, bile flow takes place peristaltically 
and the common bile duct is porous. The results presented are applicable to bile flow in the 
pathological state, where stones are formed in the bile. Evidence of slip velocity in bile flow has 
also been mentioned in Section 1. An attempt has been made to investigate to some important 
bio-fluid dynamical phenomena associated with peristaltic transport. 

Quantitative estimates have been presented for the axial velocity and critical pressure for 
reflux as a function of Reynolds number, pressure gradient, porosity parameter, Darcy num- 
ber, slip parameter, amplitude ratio as well as wave number. The results and observations are 
found to be in good agreement with those reported in [2, 33, 34, 35]. 

From the observations of this study, it reveals that the velocity proflle strongly depends on sev- 
eral parameters, viz. the mean pressure gradient, porosity, Darcy number as well as amplitude 
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ratio. The study further indicates that the reflux will occur if pressure gradient attains a certain 
critical value. 

The present theoretical investigation motivated towards the peristaltic transport of bile in 
the bile duct in the presence of stones can serve as a model that bears the potential to enrich 
our understanding of the related physiology of the problem. The important conjectures that 
can be made out of the study are as follows : 
(i) Bile velocity decreases as the stones increases. 

(ii) When bile contains a very large number of stones, reflux occurs when the critical pressure 
is quite small. 
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